Robust extended states in a topological bulk model with even spin-Chern invariant 
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This paper demonstrates the existence of topological models with gapped edge states but pro- 
tected extended bulk states against disorder. Such systems will be labeled as trivial by the current 
classification of topological insulators. Our finding calls for a re-examination of the definition of a 
topological insulator. The analysis is supported by extensive numerical data for a model of non- 
interacting electrons in the presence of strong disorder. In the clean limit, the model displays a 
topological insulating phase with spin-Chern number C s =2 and gapped edge states. In the pres- 
ence of disorder, level statistics on energy spectrum reveals regions of extended states displaying 
levitation and pair annihilation. Therefore, the extended states carry a topological invariant robust 
against disorder. By driving the Fermi level over the mobility edges, it is shown that this invariant 
is precisely the spin-Chern number. The protection mechanism for the extended state is explained. 

PACS numbers: 73.43.-f, 72.25.Hg, 73.61.Wp, 85.75.-d 



According to a widely accepted definition pQ, a topo- 
logical insulator is a material that does not conduct elec- 
tricity in the bulk but displays dissipationless conducting 
channels at the edges. This definition has been recently 
debated based on new topological models [2 , 3 that dis- 
play gapped edge states and yet they cannot be connected 
to a trivial insulator without closing the insulating gap 
or breaking the symmetry that defines them. These new 
models display anomalous responses and properties that 
are topologically protected, but they don't fit in the cur- 
rent universal classification scheme [4 j. The problem of 
classifying the insulators has been reopened. 

The present paper gives hard evidence that indeed, the 
criterium based solely on the edge states is too restrictive, 
and that large classes of materials with interesting and 
potentially useful protected properties can be neglected 
because of that. We devised a 2-dimensional (2D) lat- 
tice model that has gapped edge states but protected ex- 
tended bulk states against disorder. The model is based 
on the Kane-Mele Hamiltonian for graphene [5j |6] , 
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but we use spin operators s appropriate for spin_J par- 
ticles (thus a = ±|,±|). The notation in Eq. 



h 2 

1 is ex- 



plained in Ref. [7 . The resulting model is suiteclior the 
present analysis because a) we can compare with the orig- 
inal Kane-Mele model, b) its bands are highly entangled 
so we are not dealing with just two copies of the orig- 
inal Kane-Mele model, and c) the model remains time- 
reversal invariant. The original and the new model will 
be referenced as the spin | and | models, respectively. 

In the clean limit, the spin | model displays eight 
bands separated by an insulating gap, which closes for 




FIG. 1. The bulk (left) and ribbon (right) bands of the spin § 
model for the C s =2 case (upper) with £=0, 77=0.6 and A=0.3 
and for the C s =0 case (lower) with £=0.6, 77=0 and A=0.3. 



exceptional values of the parameters (£, 77, A) on a surface 
that separates the 3D parameter space into distinct re- 
gions, among which one with C s =2 and one with C s =0. 
For example, C s =2 for £=0, 77=0.6 and A=0.3. A plot of 
the bulk bands for this case is shown in Fig. [TJa). The 
bands for a ribbon geometry with open boundary condi- 
tions, plotted against the conserved momentum, is shown 
m Fig.gb). There are four edge bands emerging from 
the bulk spectrum, which cross each other at k points 
where the Kramer's degeneracy is not protected, so the 
bands hybridize at these crossings and become gapped. 
Note that the Rashba term responsible for this hybridiza- 
tion is large (the avoided crossings occurs far away from 
k=0) yet the gap is small. If we choose £=0.6, 77=0 and 
A=0.3, then C s =0. A plot of the bulk and ribbon bands 
are shown in Figs[TJc-d). The ribbon bands display an 
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FIG. 2. Energy spectrum and the variance of the level spac- 
ings for the spin | model. The top corresponds to the topo- 
logical case £=0, 77=0.6 and A=0.3 with C s =l and the bottom 
to the trivial case t=0.6, ?7=0 and A=0.3 with C s =0. 

almost clean bulk gap, except for accidental edge bands 
lying very close to the edges of the bulk spectrum. 
We now add spin-indepedent disorder: 

Ho -^Hu = # + W^2u n \n,a)(n,a\, (1) 
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where u n are random entries uniformly distributed be- 
tween — \ and - . The amplitudes u n are the same for all 
the states |n, a) in a unit cell. We use the level statistics 
analysis to probe the localized/delocalized character of 
the bulk quantum states, which was previously demon- 
strated to be extremely effective for Chern insulators [8] 
and the spin \ model [7 . It involves an exact diagonal- 
ization of (and later of P^&P^) on a large lattice with 
periodic boundary conditions, and a large number (10 3 ) 
of disorder configurations. Energy levels are collected 
from a small window around a given energy e and the 
level spacings between the collected levels are computed. 
By repeating this procedure for all disorder configura- 
tions, we generate an ensemble of level spacings which 
is statistically analyzed. We compute histograms show- 
ing the distribution of the level spacings and the vari- 
ance {s 2 ) — (s) 2 for various energies e. The level spacing 
distributions are compared with the appropriate Wigner- 
Dyson surmise distributions, P G vE(s) = ^-e~^ s for the 

unitary case and P G se=^3 s A e~^ s<2 for the symplec- 
tic case, and to the Poisson distribution. If there is 
an agreement with the Wigner-Dyson surmise, which re- 
mains unchanged as the lattice size is increased, then 
one can safely conclude that the states at that energy 



are delocalized [9j[10]. In such cases, the variance of the 
level spacings converges to the theoretical variance of the 
Wigner-Dyson distributions, namely, 0.178 if the opera- 
tor is in the unitary class and 0.104 if in the symplectic 
class. When the distribution of the level spacings agrees 
with the Poisson distribution and the variance is large, 
one can safely conclude that the states are localized. 

Throughout this paper, we plot the spectrum of ran- 
dom operators as follows. For each disorder configura- 
tion, we place the set of eigenvalues on a vertical line 
and then we bring these lines near each other to generate 
plots like the ones seen in Fig. [2] The representation is 
useful because one can see with naked eyes the regions 
where the eigenvalues display large fluctuations, such as 
the fuzzy regions near the edges of the spectra in Fig. [2j 
which is an indicative of localization. Sometime one can 
also see the level repulsion (see Fig. 3 of Ref. [8 ), which is 
an indicative of derealization. On top of the spectrum, 
we place the variance of the level spacings (always on a 
scale from to 1) so that we can corroborate its values 
with the fuzziness of the spectrum. 

When decifering the data, it is instructive to start 
from the spin | model, which was extensively analyzed 
in Ref. [7]. In the first panel of Fig. [2J corresponding 
to a moderate disorder iy=3, one can see energy re- 
gions where the variance is large but also two distinct 
regions where the variance becomes extremely close to 
0.104. This feature was shown [7] to remain unchanged 
when the lattice size is increased so one can safely con- 
clude that these regions contain extended states. This 
finding is in line with previous studies based on trans- 
fer matrix analysis [TT], Chern parity [12] or the non- 
commutative Z2 invariant [13], and it is widely accepted 
nowadays. As the disorder is increased, the delocalized 
spectral regions in Fig. [2] don't suddenly disappear but 
instead they migrate towards each other until they meet 
and only then they disappear. This behavior is called lev- 
itation and pair annihilation and it is the hallmark of the 
extended states carrying a topological number. Such ex- 
tended states cannot disappear unless they collide with 
other extended states carrying the opposite topological 
number, in which case the topological invariants cancel 
each other and the states become trivial and they im- 
mediately localize when the disorder is further increased. 
If we look instead at the trivial insulator in Fig. [2j the 
behavior is completely different: the levitation and pair 
annihilation is absent and the spectrum becomes local- 
ized even at small disorder. 

The data for the spin | model is presented in Fig. [3] 
The calculations for this case are much more tedious be- 
cause the dimension of the Hilbert space doubles when 
we go from spin \ to spin | . Consequently, we were able 
to complete the calculations only for 20x20, 25x25 and 
30x30 unit cells lattices. Even for these lattices, the com- 
putational effort was substantial, measuring in months 
of CPU time. Nevertheless, we saw a convergence with 




the lattice size so our conclusions are robust. The first 
panel of Fig. [3] corresponds to the topological case (whose 
clean limit was analyzed in Fig. [TJa-b)) with moderate 
disorder W=3. The histograms of the level spacing en- 
sembles recorded at various energies are shown in the 
panels immediately to the right of the energy spectrum. 
The distributions are Poisson near the edges of the spec- 
trum, while they become similar to P G se(s) as we move 
towards the center of the bands. The variance, computed 
for the three lattice sizes mentioned above, has regions 
where it becomes practically equal to 0.104, the variance 
of Pgse(s). These features are stable as the size of the 
lattice is increased, hence we can safely infer the occur- 
rence of delocalized bulk spectrum. As the disorder is 
increased, the delocalized spectral regions levitate until 
they touch and then disappear at around W=ll, totally 
analogous to what we have seen in the spin \ calculation. 
There is only one possible conclusion: for the topologi- 
cal case, the spin | model has protected extended states 
carrying a non-trivial topological invariant. If we exam- 
ine the trivial case (whose clean limit was analyzed in 
Fig. [TJc-d)), we see that the levitation and pair annihi- 
lation is absent and instead the states localize even at 
moderate disorder, such as W=5. 

We are going to show in the following that the bulk 
extended states seen for the spin | model are protected 
by the spin-Chern number introduced in Ref. [14 , us- 
ing the new formulation given in Ref. [15]. Let a be 
defined by <r|n, a) = sgn(a)|n, a). In general, the spec- 
trum of PuCrPuo is symmetric relative to the origin with 
the positive and negative parts separated by a mobility 
gap [3 [16] . If P^ define the projectors onto the pos- 
itive/negative spectrum of P u aP u , then one can define 
corresponding Chern numbers: 

C±= 2 ~f (E.M^i, P±],[* 2,i£]]|0, a)), (2) 

via the non- commutative formula of Bellissard et al [T7] . 
The outer angular parenthesis in Eq. [2] signify disorder 



average. Then the non- commutative spin-Chern number 
is defined as C S = ^[C+ — C_]. 

In Ref. [15] it was realized that the projectors P^- fit 
into the non-commutative theory of the Chern number 
[T7] and consequently we can assert with absolute con- 
fidence that C±, and therefore C s , are quantized and 
invariant as long as Pj 1 remain localized. The derealiza- 
tion of these projectors can happen via two mechanisms: 
1) the mobility gap of P^aP^ closes, or 2) the projector 
Pu itself becomes delocalized. For the spin \ model, the 
mobility gap of P^aP^ was found to be extremely robust 
and only the second mechanism was observed [7] . We are 
going to examine the numerical data and use analytic ar- 
guments to show that same happens for the | model. If 
that is the case, then a nonzero value of C s would nec- 
essarily imply the existence of extended bulk states. In- 
deed, imagine that we continuously lower the Fermi level 
Ep until it reaches the bottom of the spectrum, where 
Puj=0 and consequently C s =0. We can then see that C s 
changes its value during this process, so for some values 
of Ep the P^ must became delocalized. Since the first 
mechanism is absent, this derealization can come only 
from the derealization of the whole projector P w , which 
implies the existence of bulk extended states. 

In Fig. [4] we mapped C s and the spectrum of P^aP^ 
as the Fermi level was continuously lowered. We com- 
puted C s by evaluating Eq. |on a finite lattice with 
periodic boundary conditions by following the fast con- 
verging procedure introduced in Ref. [8] and detailed in 
Ref. [7 . We mention that none of the calculations pre- 
sented here would have been possible without this proce- 
dure, which increases the accuracy and reduced the CPU 
time tremendously when compared with the traditional 
twisted boundary conditions calculations [12 . Let us 
concentrate on the spectrum of P^aP^ first. This opera- 
tor is in the unitary class, so if mechanism 1) would have 
taken place, one should see in the spectrum of P^cfP^ a 
very narrow band of extended states moving towards the 
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FIG. 4. Bottom: The spectrum of H u and the variance of the energy level spacings at W—?>. Middle: The spectrum of P w a"P w 
and the variance of the level spacings for various Fermi levels. Upper: The numerical value of the spin-Chern number for the 
marked Fermi levels. The calculation was completed on a 30x30 lattice. 



origin and closing the mobility gap. Such narrow bands 
of extended states would have had a distinct signature 
on the level statistics analysis [8] and the variance would 
have had sharp valleys where it abruptly converged to 
0.178. These valleys should move towards the origin in 
order to close the mobility gap of P^aP^. Clearly we 
don't see that in the data of Fig. [4j Instead, we see a 
sudden total collapse of the variance onto the value 0.178. 
This collapse happens immediately after the Fermi level 
enters the energy region where the variance of the energy 
level spacings becomes 0.104. Such collapsing behavior 
was previously seen in the entanglement spectrum of a 
Chern insulator [8] when the Fermi level crossed the An- 
derson transition point. The origin of that collapse is 
well understood now and comes from the derealization 
of P u . Similarly for the present data, the collapsing be- 
havior can only come from the derealization of P u . 

The data on the spin-Chern number strongly supports 
the above conclusion. As we lower the Fermi level, C s 
is seen to remain quantized at 2 until the Fermi level 
reaches the point where the variance of the energy level 
spacings becomes equal to 0.104, and where we infer that 
the protected extended spectrum starts. From there on, 
C s starts decreasing until Ep exists the energy region 
of extended states, when C s sets to exactly the value 0. 
Putting all three data analysis together, the energy level 
statistics, the P^aP^ level statistics and the C s compu- 
tation we have no alternative but to conclude that the 
spin | model has protected extended bulk states and the 
protection is provided by the spin-Chern invariant. 

One question remains, why does the mobility gap of 
PujCtPuj remain open? There are two reasons for this. 
First, the disorder is spin independent, so it has negligi- 



ble direct effects on P^aP^. It is reasonably clear that 
the mobility gap of the Hamiltonian is more sensitive to 
disorder than the mobility gap of P^cfP^ is. Second, the 
following identity: 

2P± =P^± P u aP u (l - (*[a,P w ]) 2 p (3) 

shows that Pj 1 relate to the Green's function of the self- 
adjoint operator i[a, P^], evaluated at ±1. The spectrum 
of i[cr, PJ\ is always inside the interval [—1, 1] so we are 
basically probing the edges of its spectrum. But as long 
as P u is localized, i[a, P u ] is a random matrix with well 
behaved entries so the edges of its spectrum are expected 
to be localized. These two arguments show that mobility 
gap of P^aP^ is extremely robust. 

In conclusion, we have presented three independent nu- 
merical arguments all converging to one conclusion, that 
of existence of protected extended bulk states in a sys- 
tem with trivial invariant. The protection mecha- 
nism steams from the insensitivity of the mobility gap of 
PuCrPu to the spin-independent disorder, fact that gives 
the spin-Chern invariant a status similar to that of the 
Chern invariant in Chern insulators. This protection is 
probably not universal, but it is clearly taking place in 
the particular model examined here and will very likely 
occur in many other models and real materials. In a 
good sense, our analysis shows how little we understand 
the disordered topological insulators and how much is 
there to explore. 
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